Entropy of entanglement and multifractal exponents for random states 
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We relate the entropy of entanglement of ensembles of random vectors to their generalized fractal 
dimensions. Expanding the von Neumann entropy around its maximum we show that the first order 
only depends on the participation ratio, while higher orders involve other multifractal exponents. 
These results can be applied to entanglement behavior near the Anderson transition. 
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Entanglement is an important characteristics of quan- 
tum systems, which has been much studied in the past 
few years due to its relevance to quantum information 
and computation. It is a feature that is absent from clas- 
sical information processing, and a crucial ingredient in 
many quantum protocols. In the field of quantum com- 
puting, it has been shown that a process involving pure 
states with small enough entanglement can always be 
simulated efficiently classically [ij. Thus a quantum al- 
gorithm exponentially faster than classical ones requires 
a minimal amount of entanglement (at least for pure 
states). Conversely, it is possible to take advantage of 
the weak entanglement in certain quantum many-body 
systems to devise efficient classical algorithms to simu- 
late them [2|. All these reasons make it important to 
estimate the amount of entanglement present in different 
types of physical systems, and relate it to other proper- 
ties of the system. However, in many cases the features 
specific to a system obscure its generic behavior. One 
way to circumvent this problem and to extract generic 
properties is to construct ensembles of systems which af- 
ter averaging over random realizations can give analytic 
formulas. Such an approach has proven successful, e.g. in 
the quantum chaos field, where Random Matrix Theory 
(RMT) can describe many properties of complex quan- 
tum systems. 

One of the interesting questions which have been ad- 
dressed in many studies (see e.g. [3| and references 
therein) is the behavior of entanglement near phase tran- 
sitions. It has been shown that the entanglement of the 
ground state changes close to phase transitions. For ex- 
ample, in the XXZ and XY spin chain models, the en- 
tanglement between a block of spins and the rest of the 
system diverges logarithmically with the block size at the 
transition point '4||, making classical simulations harder. 
However, such results cannot be applied directly to sys- 
tems where the transition concerns one-particle states, 
for which entanglement has to be suitably defined. A 
famous example is the Anderson transition of electrons 
in a disordered potential, which separates localized from 
extended states, with multifractal states at the transition 
point. Previous works [5| have described the lattice on 
which the particle evolves as a spin chain and studied en- 
tanglement in this framework. However, the lattice can 



alternatively be described in terms of quantum computa- 
tion with a much smaller number of two-level systems [6| . 

In this paper, we study entanglement of random vec- 
tors which can be localized, extended or multifractal in 
Hilbert space. We consider entanglement between blocks 
of qubits. In the case of the Anderson transition, this 
amounts to directly relate entanglement to the quantum 
simulation of the system on a rir-qubit system, the num- 
ber of lattice sites being 2"'' rather than n^ as in [5|]. 
Entanglement of random pure states was mainly stud- 
ied in the case of columns of matrices drawn from the 
Circular Unitary Ensemble (CUE) [7|]. However, such 
vectors are extended and cannot describe systems with 
various amounts of localization, from genuine localiza- 
tion to multifractality. Recently it was shown in [8|, |9| 
that for localized random vectors, the linear entangle- 
ment entropy (first order of the von Neumann entropy) 
of one qubit with all the others can be related to the 
localization properties. Here we develop this approach 
to obtain a general description of bipartite entanglement 
in terms of certain global properties for random vectors 
both extended and localized. First we show that for any 
bipartition, the linear entropy can be written in terms 
of the participation ratio, a measure of localization. We 
then show that higher-order terms also depend on higher 
moments of the wavefunction. In particular, for multi- 
fractal systems they are controlled by the multifractal 
exponents. 

Bipartite entanglement of a pure state {ip) belonging to 
a Hilbert space Ha^^T^b is measured through the entropy 
of entanglement, which has been shown to be a unique 
entanglement measure jlO| . Let pA be the density matrix 
obtained by tracing subsystem B out of p = |V')(V'|- The 
entropy of entanglement of the state with respect to the 
bipartition {A, B) is the von Neumann entropy of pA, 
that is S* = — tr (pa ^og2 Pa)- It is convenient to define 
the linear entropy as Sl = ^jri^-'^ ~ trp^), where d — 
diraTiA < dimTis- The scaling factor ensures that Sl 
varies in [0,1]. We will show that the average value of 
S over a set of random states can be expressed only in 
terms of averages of the moments of the wavefunction 
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provided some natural assumptions are made. Here we 
consider ensembles of random vectors of size N = 2'^'' 
with the following two properties: i) the phases of 
the vector components are independent, uniformly dis- 
tributed random variables, and ii) the joint distribution 
P{xi, . . . ,xn) oi the modulus squared of the vector com- 
ponents is such that all marginal distributions P{xi), 
P{xi,Xj) for i y^ j, P{xi,Xj,Xk) ioT i ^ j ^ k and so 
on, do not depend on the indices. As a consequence, 
all correlators (IV'ii P*MV'i2p^^ • • ■) of the components of 
\ip) are independent of the indices i\ ^ 12 ^ ■ ■ ■ involved. 
Random vectors realized as columns of CUE matrices are 
instances of vectors having such properties. 

Let us first consider the simplest case of entanglement 
of one qubit with respect to the others. Then d — 2 
and the linear entropy Sl is simply the tangle r, or the 
square of the generalized concurrence [ll| . It is given by 
T — Adet PA- If we consider a vector \ij}) of size N, the 
bipartition with respect to qubit i splits the components 
ijjj of IV') into two sets, according to the value of the ith 
bit of the binary decomposition of j. If IV''-^-') and IV-''-^-') 
are the two corresponding vectors, the linear entropy is 

T = 4 ((^(o)|V'(o))(^W|V'(i)) - KV'^^^IV'^'^)!') . (2) 

After averaging r over random phases, only the diagonal 
terms survive in the scalar product \{^^^'>\ip'-^'>)\'^. Since 
it is assumed that two-point correlators of the vector j?/') 
do not depend on indices, their average can be expressed 
solely in terms of the mean moments, as (IV'iPlV'iP) = 
N(N-i) f'-"' * 7^ J- Normalization of \ip) implies pi = 1. 
As vectors jV''-'''') and \ij^^^) always contain components 
of \ip) with different indices, we get 
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Since p2 — 1/^ where ^ is the inverse participation ratio 
(IPR), Eq. © is exactly the Eq. (3) of Ref. [9]. Let us 
now turn to the general case, and consider the entropy of 
entanglement of v qubits with n^ — v others {{v, Ur — v) 
bipartition). The vector 1-0) is now split into vectors 
IV'^-'''), < J < 2"^ — 1, depending on the values of the v 
qubits. The reduced density matrix pA then appears as 
the Gram matrix of the IV" )j and the linear entropy is 



2"^- 1 



2"-! 

E 



|(V'«|0(^))|2 



(4) 



When averaging over random vectors, each term in 
Eq. (|4]) with i ^ j yields 2""'~'' two-point correlators, 
while each term with i = j yields 2""'~'^(2"''~'^ — 1) two- 
point correlators and 2"''^'' terms of the form (|V'i|''). 
Inserting these expressions into ([!]) gives {Sl) = {N — 
2'')(1 - {P2))/{N - 1), which generalizes Eq. 0. The 



first-order series expansion of the mean von Neumann 
entropy around its maximum can be expressed as 



(S) 



1 



21n2 



1 



N -2" 

N- 1 



(1 - (P2)) 



(5) 



with p2 = 1/S,- Equation ([5]) shows that for any par- 
tition of the system into two subsystems, the average 
bipartite entanglement of random states only depends at 
first order on the localization properties of the states, 
through the mean participation ratio. For CUE vectors, 
formula (O reduces to the expression for the mean en- 
tanglement derived earlier in [13]. More interestingly, 
this formula also applies to multifractal quantum states. 
There, the asymptotic behavior of the IPR is governed 
by the fractal exponent D^, where one defines generalized 
fractal dimensions Dq through the scaling of the moments 
Pq oc N~^''^'i~^\ Thus the linear entropy is only sensitive 
to a single fractal dimension. These results imply that 
entanglement grows more slowly with the system size for 
multifractal systems. 

To test the relevance of Eq. ([5|) for describing entan- 
glement in realistic settings, we consider eigenvectors of 
N X N unitary matrices of the form 
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where (j)k are independent random variables uniformly 
distributed in [0, 27r[. These random matrices display in- 
termediate statistical properti es Il3l . and possess eigen- 
vectors that are multifractal [iJl, both features being 
tuned through the value of the real parameter 7. We also 
illustrate Eq. ^ with eigenstates of a many-body Hamil- 
tonian with disorder and interaction H = 'Y^^Tia^ + 
J2i<i Jij^i^] ■ This system can describe a quantum com- 
puter in presence of static disorder [l5| . Here the ai are 
the Pauli matrices for qubit i, energy spacing between 
the two states of qubit i is given by F^ randomly and uni- 
formly distributed in the interval [Aq— (5/2, Aq -1-5/2], and 
the Jij uniformly distributed in the interval J— J, J] repre- 
sent a random static interaction. For large J and 5 « Aq, 
eigenstates are delocalized in the basis of register states, 
but without multifractality. They display properties of 
quantum chaos, with eigenvalues statistics close to the 
ones of RMT [15|. For both systems (unitary matrices 
and many-body Hamiltonian) , components of the eigen- 
vectors have been shuffled in order to reduce correlations, 
but leaving the peculiarities of the distribution itself un- 
altered. Figure [1] plots the first-order expansion ([5]) as 
a function of the mean IPR for three different biparti- 
tions, showing remarkable agreement with the exact {S), 
both for multifractal (Fig. [U left panel) and non fractal 
(right) states, and even for moderately entangled states. 
The agreement is better for the non-fractal system than 
for the multifractal one. This can be understood from 
the study of higher order terms in the entropy. 




10 20 30 40 50 CO 



120 160 



(1/0-^ 



(1/0-^ 



FIG. 1; (Color online) Mean entropy of entanglement as a 
function of the mean IPR. Left: eigenvectors of (|6]) with 
7 = 1/3; the average is taken over 10^ eigenvectors. Right: 
eigenvectors of the Hamiltonian H (see text) with 5 — Ao and 
J/5 = 1.5; average over N/16 central eigenstates, with a total 
number of vectors « 3 x 10^. Triangles correspond to «/ = 1, 
squares to ;/ = 2 and circles to i/ = nr/2, with n^ = 4 — 10 
(bipartition of the i/ first qubits with the nr—u others). Black 
symbols are the theoretical predictions for (S) at first order 
(Eq. (O) and green (grey) symbols are the computed mean 
values of the exact (S). 



Indeed, while the hnear entropy does not depend on 
other fractal dimensions than D2, the entropy of entan- 
glement does. If we go back to the case of a (l,nr — 1) 
bipartition of the system, the entropy of entanglement 
can be expressed in a simple way as a function of t as 



S{t) = h 
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where h(x) = — a;log2 x — {1 — x) log2(l — x). The series 
expansion of S{t) up to order ?7i in (1 — r) reads 
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The tangle r corresponds, up to a linear transformation, 
to Si{t). Let us now calculate the average of higher 
orders in this expansion. The second-order expansion of 
S{t) involves calculating the mean value of 
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where the star denotes complex conjugation and Ui, Vi are 
the components of \il)^'^'^), IV'*'^'') respectively. Under the 
assumption of random phases, only terms whose phases 
cancel survive in Q . Since the phases of all components 
of \^) are independent, cancelation of the phase can only 
occur if the sets {i, k} and {j, 1} are equal. Thus 



The correlators in Eq. pO|) can be expressed as a func- 
tion of the moments as follows. Using standard notations 
[la |. we will denote by A h n a partition A = (Ai, A2, . . .) 
of n, with Ai > A2 > .... For any partition A h n, we 
define px — PA1PA2 ■ • -i where px. are given by (H]). The 
monomial symmetric polynomials are defined as max = 
X^dV'iP'^^ |'/'2p^^ • ■ •) (the sum runs over all Vx permuta- 
tions of the Ai), and we set cx = rnx/Vx- The px and mx 
are related by the simple linear relation pA = X)u Lx^im^, 
where Lx^ is an invertible integer lower-triangular matrix 
([16j, p. 103). Upon our assumption ii), any correlator of 
the form (|i/'ii P''MV'J2p''^ • ■ ■) is equal to a cx for some 
partition A of n, and thus can be expressed as a func- 
tion of the moments. For instance the two correlators in 
Eq. (fTU]) are respectively equal to Cim and C22. Treating 
similarly all terms involved in t^ gives 

(r^) = iV(iV-2)(iV2-6iV+16)ciiii (11) 

+ 4iV(iV-2)(7V-4)c2ii+4Ar(Ar-2)c22. 

This term involves the calculation of three correlators. 
Using the relation between the cx and px and the fact 
that the vectors are normalized to one we get 

{pD - (P4) {P2) - {pD - 2(P3) + 2(P4) 

C22 = -TTnTr TV' C211 = 



Clin = 



N{N-l)' N{N ~1){N -2) 

1 - 6(p2) + 8(p3) + i{pl) ~ GjPi) 
N{N-l){N-2){N ~3) 



(12) 



The calculation of the general term (r") can be per- 
formed along the same lines. Expanding ^ we get 



r" = 4"^ ("j(_l)"-fc 




(13) 
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The expansion of I '^- u*ViUjV* ) contains products of 

the form (u^^ . . .Ui^)*Uj-^ . . .Uj^. Only terms where the 
phases coming from the u* compensate those coming 
from the Uj^, survive when averaging over random phases. 
Thus we keep only terms where {ji, . . . , jt} is a permuta- 
tion of {«!, . . . , i(}. If Vk is the number of permutations 
of a set K , the average of (fT3l) over random vectors reads 
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Terms with the same correlator can be grouped together. 
Each correlator in P^ is some caua' , with A, A' parti- 
tions of n. For A h n and ^ \~ k we define the coeffi- 
cient Axu = TT y^ /"~. Si' , where the sum runs over all 
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FIG. 2; (Color online) Relative difference of the entropy of 
entanglement ((Tjl and its successive approximations Sm {m = 
1, 2) with respect to the number of qubits for eigenvectors of 
((6]) for (left) 7 = 1/3 and (right) j = 1/7. The average is 
taken over lO'' eigenvectors, yielding an accuracy < 10~® on 
the computed mean values. Green triangles correspond to the 
first-order expansion Si, blue squares and red circles to the 
second-order expansion 52. The difference between the latter 
two is that for blue squares (pi) appearing in Eq. (|12|) has 
been replaced by (p2)^ yielding a less accurate approximation. 
Dashed line is a linear fit yielding 1 - (Si)/(S) ~ JV""*** for 
7 = 1/3 and N'^-^^ for 7 = 1/7. 



vectors s = (si, . . . , Sjv/2) which are permutations of A, 
and k — {p,i, ... ,^1^/2). We have used the notations 
a ~ (oi, a2, . . .) and a}. — ai\a2\ ■ . .. Finally we get 
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(15) 
which provides an expression for the nth order for (5) 
as a function of the (pa)- In the special case of CUE 
random vectors, by resumming the whole series we re- 



cover after some algebra the well-known result [17 1 



(5(r)) = j^ X]fc^/2+i \- Note that similar expressions 
can be derived for a general (y, Ur — v) bipartition. In this 
case, the entropy S = — tr(p^log2 pa) can be expanded 
around the maximally mixed state po = 1/2"^, as 
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After averaging over random vectors, one can check that 



the traces in (fTB)l can be written as {tip'} 



l^Xhk "a 



with a\^ some integer combinatorial coefficient. The en- 
tropy can thus be written as a linear combination of c\ 
with rational coefficients that can be expressed in terms 
of the a\^' . 

In Fig. [5] we illustrate the accuracy of higher-order 
terms in the series expansion of S for multifractal random 
vectors by comparing the first and second-order expan- 
sion for eigenvectors of the matrices ([6]). As expected, the 
second-order expansion is much more accurate than the 
first order one and gives a much better estimate of the 
mean entropy of entanglement already for small system 



sizes. For large N, the dominant term in S2 is ex (p^)- 
Numerically we obtained (^2) ^ N~'^-^^ for 7 = 1/3, and 
(pf) ^ Af^i-53 fQj- ry — 2/7, which is indeed consistent 
with the slopes of the finear fit of log2(l — {Si)/{S)) (see 
Fig. [2]). If one replaces (p|) appearing in Eq. p^ by 
(p2)^ (squares in Fig. [5]), the second-order expansion is 
now governed only by three multifractal dimensions D2, 
D3, D4. Although it becomes less and less accurate with 
the system size because of the increase of the variance of 
P2 , it remains a very good improvement over the first or- 
der in the case of moderate multifractality (Fig. [21 right). 

Our results show that the entanglement of random vec- 
tors directly depends on whether they are localized, mul- 
tifractal or extended. The numerical simulations for dif- 
ferent physical examples show that our theory describes 
well individual systems whose correlations are averaged 
out. Previous results I9| have shown that Anderson- 
localized states have entanglement going to zero for large 
system size. The present work shows that multifractal 
states, such as those appearing at the Anderson transi- 
tion, approach the maximal value of entanglement in a 
way controlled by the multifractal exponents. Although 
extended and multifractal states are both close to maxi- 
mal entanglement, the way multifractal states approach 
the maximal value for large system size is slower. 
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